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ABSTRACT 



Aims. We analyze 19 bright variable X-ray sources detected in the XMM-Newton Extended Survey of the Taurus Molecular Cloud 
(XEST), in order to characterize the variations with time of their coronal properties and to derive informations on the X-ray emitting 
structures. 

Methods. We performed time-resolved spectroscopy of the EPIC PN and MOS spectra of the XEST sources, using a model with one 
or two thermal components, and we used the time evolution of the temperatures and emission measures during the decay phase of 
flares to derive the size of the flaring loops. 

Results. The light curves of the selected sources show different types of variability: flares, long-lasting decay or rise through the whole 
observation, slow modulation or complex flare-like variability. Spectral analysis shows typical quiescent plasma temperatures of ~ 5 - 
10 MK and ~ 15-35 MK; the cool component generally remains constant, while the observed flux changes are due to variations of 



the hot component. During flares the plasma reaches temperatures up to 100 MK and luminosities up to 
inferred from flare analysis are generally smaller than or comparable to the stellar radius. 
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1. Introduction 

Young pre-main sequence (PMS) stars are well known to exhibit 
. £h ' stxon g X-ray emission at levels orders of magnitude higher than 
K> , generally observed in older active stars (|Feigelson & Montmerlel 
Q ■ ll999tlFavata & Micelall2003fc lGudelll2004l) Furthermore. X-rav 
' time variability is almost always detected in PMS stars, very of- 
ten in the form of flares, with a shape similar to solar flares, or 
of slow modulation due e.g. to rotation (Feigelson & Mont merlel 
1999; IWolk et alj|2005 IFlaccomio et alj|2005l) . On the Sun, the 
X-ray time variability is due to flare activity on timescales of a 
few hours, to the appeareance of new active regions combined 
with the surface rotation (timescales of a few weeks) and to the 
11-yr solar activity cycle. The flare activity observed in PMS 
stars suggests that a scaled version of the solar corona should 
be present in these objects. The plasma magnetically confined in 
loop-like structures would be heated by mechanisms similar to 
those present in the solar corona, and would release the energy 
in the form of flares. However, observations show that Classical 
T Tauri stars (CTTS), characterized by accretion from the cir- 
cumstellar disk, are statistically X-ray under-luminous with re- 
spect to Wea k-lined T Tauri stars (WTTS), where disks are ab- 
sent or weak dStelzer & Neuh auser 200 it IFlaccomio et al.f 2003; 



Preib isch et ail 120051: iFranciosini etafl 120061: iTelleschi et all 
2006bl) . M oreover, high-resolut i on spectroscopy of the C TTS 
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TW Hv a dKastner et al.1 l2002t IStelzer & Schmitll |2004|) and 
BP Tau ( Schm itt et alj|2005l) indicate the presence of very high 
densities (n e ^ 10 12 cm 4 ) and a dominant cool plasma com- 
ponent at 3 MK in TW Hya, that have been attributed to emis- 
sion from plasma heated by an accretion shock. Enhanced X- 
ray emission and/or significant spectral variations have been ob- 
served during accretion outbursts in VI 647 Ori and VI 1 18 Ori 
(Kast ner et al] 12004 120061: iGrosso et all l2004t lAudard et al.1 
120051) . These results suggest that accretion may play a role in 
the X-ray emission process, either by influencing the magnetic 
structure, or by providing an alternative X-ray production mech- 



The study of the time variability of the X-ray emission from 
PMS stars allows us to gain insights into the structure and the 
heating mechanisms of stellar coronae. In particular, the analysis 
of flares constitutes a diagnostic t ool to infer the size of the X- 
ray emitting structures. Recently, F avata etail d2005l) studied a 
sample of intense flares observed in the Orion Nebula Cluster 
as part of the Chandra Orion Ultradeep Project (COUP), finding 
that flares on PMS stars occur both in small loops, of size less 
than a stellar radius, similar to older active stars, and in large 
loops, up to 10 - 20 R+, likely connecting the stellar surface and 
the circumstellar disk. 

In this paper we study the X-ray properties of a sample of 
PMS stars of the Taurus Molecular Cloud (TMC) showing sig- 
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nificant time variability. The TMC is one of the nearest regions 
of star formation (d = 140 pc), containing ~ 340 known mem- 
bers in an area ~ 100 square degrees large, and is characterized 
by a low stellar density (1 - 10 pc -3 ), by the lack of massive stars 
and by a significant fraction of binary or multiple systems. Star 
formation appears to have occurred in several epochs during the 
lastlOMyr. 

Our work complements the study by IStelzer et al.l (|2006a), 
who analyze the X-ray time variability of TMC sources from 
a statistical point of view, and derive the frequency and energy 
distribution of flaring events. In our paper we concentrate on the 
time-dependent spectral analysis of the brightest sources, show- 
ing both flares and other kinds of variability, in order to inves- 
tigate the changes of the plasma parameters and to derive infor- 
mation on the X-ray emitting structures. 

The paper is organized as follows. In Sect. [2] we describe 
the target selection and the data analysis. In Sect. [3]we present 
the results obtained for TMC members, whereas in Sect. |4] we 
present the variability of three sources not related to the TMC. 
Discussion and conclusions are given in Sect. [5] 



2. Source selection and analysis 

2.1. Observations 

Our study is based on the data obtained from the XMM-Newton 
Extended Survey of the Taurus Molecular Cloud (XEST), a wide 
X-ray survey performed with the XMM-Newton satellite aimed 
at studying the pro perties of young PMS stars in the TMC 
dGiidel et a l. 2006a). The survey consists of 19 fields of ~ 33 ks 
exposure each, covering the densest regions of the cloud, plus 9 
archival exposures lasting up to ~ 130 ks; the total area covered 
by the survey is about 5 square degrees. 

A detailed description of the survey, including details on the 
primary reduction of the raw EPIC MOS and PN datasets, on the 
source detection procedure and on the definition of the extraction 
regions for t he source and backg round light curves and spectra 
is reported in lGiidel et alj d2006al) . 

The work presented here uses data from all XEST fields 
except for X EST-01 and XEST- 25 that are dedicated to sepa- 
rate projects (Giide l et alj|2006bl Grosso et al., in preparation). 
We also excluded the classical T Tauri star SU Aur (source 
XEST-26-067) that is d iscussed in detail in a separate paper 
dFranciosini et alj [20 07). We also note that two of the fields, 
namely XEST-23 and XEST-24, correspond to two consecutive 
observations with the same pointing position, with a duration of 
~ 70 and 40 ks, respectively, separated by ~ 5 ks. Sources de- 
tected in both exposures have been treated as a single source in 
the following analysis, with the two exposures added together. 

2.2. Analysis of light curves and sample selection 

To select the sample of sources for this study, we first extracted, 
for all detected X-ray sources, the PN and MOS source and back- 
ground photons in the ene rgy band 0.3 -7.3 k eV, using the ex- 
traction regions defined bv lGiidel et al.l d2006al) . In order to have 
a continuous time coverage for each observation, without gaps 
in the light curves, we chose not to filter out the periods of high 
background count rate due to proton flares. We then limited our 
analysis to sources with at least 1500 net counts in the PN expo- 
sures (or in the MOS exposures in the cases where the PN was 
not available), in order to have sufficiently high count statistics 
for good time-resolved spectral analysis. This selection neces- 



sarily biases our study towards the X-ray brightest and likely 
most active sources. 

For these sources, we extracted the PN (or MOS) light 
curves, and applied to them the Maximum Likeli hood Blocks 
(MLB) method described by IStelzer et all d2006ah . in order to 
identify variable sources and to define the intervals to be used 
for time-dependent spectroscopy. The MLB algorithm divides 
the full exposure into time intervals (blocks) where the source 
is assumed to be constant; the boundaries of the time blocks are 
chosen such that between two consecutive blocks the mean count 
rate changes by more than a given significance threshold. The 
significance threshold for the change points, at the 99% level, 
was determined through simulations of constant rate light curves 
(Flaccomio et al., in preparation). We set a minimum number of 
750 net source photons for each block, in order to have enough 
counts in each spectrum to perform a reliable spectral analysis. 
We note that this choice does not allow us to examine rapid vari- 
ations, small-amplitude flaring events and low-level variability. 

The variability of the background in many XEST fields was 
significant and was taken into account using the following proce- 
dure: the MLB algorithm was used to split the background light 
curve into blocks of constant count rate level; for each of these 
blocks, the background level was scaled to the source extrac- 
tion area to yield the number of expected background photons in 
the source area for the given time interval. The resulting num- 
ber of photons is then removed from the source event file uni- 
formly across each background block. The result is a "source- 
only" event file, i.e. a background-subtracted photon time series 
for the source. 

Using the MLB results, we restricted our sample by consid- 
ering only those sources whose light curves showed significant 
variability, i.e. those that have been divided by the algorithm in 
at least two blocks with mean count rates differing by more than 
3cr. This selection process led to a sample of 16 variable X-ray 
sources associated with known PMS members of the TMC, plus 
additional three sources not related to the region or that are not 
confirmed members of the TMC. The weakest source in the final 
sample is JH 108, with ~ 2500 net PN counts. The main proper- 
ties of the selected sources are given in Table Q] The three non- 
members will be discussed separately in Sect. [4] 

The 16 TMC members are mostly late-K and early-M stars; 
five of them are classified as CTTS, ten as WTTS, and one is a 
Herbig Ae star. Six of these sources (including the Herbig Ae 
star) are binary or multiple systems unresolved in X-rays with 
XMM-Newton. 

2.3. Time-resolved spectroscopy 

We have performed a time-dependent spectral analysis of the se- 
lected sources, in order to study the changes with time of the 
plasma characteristics, i.e. temperature and emission measure 
(EM), that give us information on the origin of the emission, 
and to derive, in the case of flares, information on the size of the 
flaring region from flare modeling. 

PN and MOS spectra for each source and each time block 
have been extracted from the source and background event files. 
For the response matrices, we used the appropriate canned re- 
sponse matrix files for PN and MOS, and the ancillary response 
files produced for each individual source (see lGiidel et alj2 006a, 
for details). The PN and MOS spectra have been rebinned in or- 
der to have at least 15 counts per bin, and have been jointly fitted 
in XSPEC v. 11.3.2, using a thermal VAPEC model with one or 
two temperature components, plus a common photoelectric ab- 
sorption component. The hydrogen column density Nh was left 
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Table 1. Properties of the sample of XEST variable sources analyzed in this paper. Stellar parameters are taken from lGiidel et al.l 
(2006a). In the last column we give a classification of the observed variability. 
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as a free parameter, except in a few cases where it was not con- 
strained by the fit, and it was therefore fixed to a value derived 
for other blocks of the same source. Abundances were kept fixed, 
with their values following a pattern derived from e stimates of 
coronal abundances of X-ray active young sta rjj (Scelsi et al.l 
120051; iTelleschi et alj2005llArgiroffi et al.l2004l) . The best-fit re- 
sults are reported in Table |2j errors are 90% confidence intervals 
for one interesting parameter. 

Figs. Q]|9] show the light curves and spectral parameters de- 
rived for TMC sources; the results for TMC non-members are 
shown in Figs. fTOlfTTl For each source, we plot in the top panel 
the PN (or MOS) background-subtracted light curve and the 
background light curve, while in the other panels we show the 
evolution of the best-fit temperatures and EMs and of the hydro- 
gen column density. The time intervals derived from the MLB 
method and used for the analysis are indicated by dotted vertical 
lines. In Figs. l2l [3l 131 and [TTI we also show, in the right panels, 
the evolution of the plasma temperature and emission measure 
used for flare analysis. 



2.4. Flare modeling 

Flare light curves contain information on heating and cooling 
parameters and indirectly on the density and geometry of the 
flaring sources. As stellar flare observations are not spatially re- 
solved, various techniques have been developed to use the light 
curves and measurable quantities such as the emission mea- 
sure or the electron t emperature t o deduce flare source param- 
eters (see review bv iGiidell [20041. Each of these methods in- 
troduces simplifications and makes various assumptions on the 
heating and cooling processes, often adopted by analogy with 
solar flares. A parti cularly stra i ghtfor ward physical model has 
been developed by Real e et all Q997), that derives the heating 
fraction and the magnetic loop length from hydrodynamic con- 



1 The abundan ce s use d, relative to the solar abundances by 
lAnders & Grevessd dl989h . are: C = 0.45, N = 0.788, = 0.426, 
Ne = 0.832, Mg = 0.263, Al = 0.5, Si = 0.309, S = 0.417, Ar = 0.55, 
Ca = 0. 195, Fe = 0. 195, Ni = 0. 195. 



siderations, using only the run of flare temperature and emis- 
sion measure. Its advantage is its simple application to stellar 
observations; it has also been tested on moderate single-loop so- 
lar flares. It is, on the other hand, not applicable to flares oc- 
curring in magnetic arcades as in gradual solar flares, where a 
large number of magnetic structures is ignited and cooling in se- 
quence. Then, both the rise and the decay time are essentially de- 
termined by the history of the heating energy release rather tha n 
by cooling physics (e.g. lPneumanlll98"2tlKopp & Polettolll98 
Co mparative studies h ave been presented in Gud el et all (1200- 
and iReale et all d2004l) for an exceptionally strong flare on the 
nearby active star Proxima Centauri. We are not in a position to 
characterize flares on T Tau star s here. We henceforth will adopt 
the model of IReale et al.l d 1997b to discuss possible systematics 
in our results. 

In the case of a sing le flaring loop dec aying freely after an 
initial impulsive heating, Seri o et all dl99ll) showed that the loop 
semi-length L is related to the thermodynamic decay time t± by 
the relation: 



3.7 x 10- 4 



(1) 



where T max is the maximum temperature of the flaring plasma. 
However, if significant heating is present also during the decay 
phase, the flare decay time will be longer than r t h, and the above 
equation would lead to an overstimate of the loop length. Using 
hydrodynamical simulations of flari ng loops taking in to account 
the presence of prolonged heating, IReale et all (1 19971) showed 
that the slope f of the flare decay path in a diagram of log T 
vs log n e (or equivalently log V-EM) is a diagnostic of the pres- 
ence of residual heating during the flare decay. The loop semi- 
length L can then be derived from £ and the observed flare decay 
timescale t\ c using the following relation: 



L = 



3.7 x 10- 4 F(0 ' 



(2) 



where F(£) is the ratio between the observed and thermody- 
namic decay times, which depends on the amount of residual 
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V410 X7 = XEST-23-050+24-042 - Bin = 1200 s 



V819 Tou = XEST-23-074+24-061 - Bin = 1200 s 
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Fig. 1. Light curves and spectral fitting results for the flaring TMC sources V410 X7 (XEST-23-050/24-042) and V819 Tau (XEST- 
23-074/24-061). In the top panels, the background-subtracted PN light curve of each source is plotted with a thick solid line, while 
the thin solid line shows the background light curve, scaled to the source area. The vertical dotted lines mark the time intervals, 
derived from the MLB analysis, used for the time-resolved spectroscopy. In the other panels we show, from top to bottom, the time 
evolution of the best-fit temperatures, emission measures and hydrogen column density. The temperatures and EMs of the cool and 
hot component are indicated by diamonds and asterisks, respectively. 



heating. T max , as in Eq. ([T]), is the maximum temperature in the 
flaring loop, that can be derived from the flare peak temperature 
Jobs obtained from the spectral fits. The expression of F{() and 
the relationship between T max and r o b s depend on the instrumen- 
tal re sponse. In the case of EPIC PN we have ( iGiardino et alj 
2006): 



= 0.13 T 



1.16 

obs 



and 



F{0 = 1.36 + 



0.51 
<T-0.35 



(3) 



(4) 



The expression of F(0 is valid only for 0.35 < £ < 1.6, 
where the lower limit corresponds to the case where the heat- 
ing timescale is comparable to the flare decay timescale, while 
the upper limit corresponds to a freely-decaying loop with no 
residual heating after the impulsive energy release. 

The above formulae are calibrated on the EPIC PN response, 
however they give a correct order of magnitude estimate of the 
loop length also in the case of MOS, given the similarity of 
the instrumental responses and the wide spectral band used. We 
have therefore applied them also to the analysis of the flare of 
HBC 427, for which PN is not available. 

3. Sources associated with TMC members 

We have divided the sources associated with TMC members into 
three classes according to the type of variability observed in 
their light curves: (i) sources showing evident flares, either with 
a "classical" shape, with fast rise and slower decay, or with a 



gradual rise and/or a symmetrical shape; (ii) sources with slow 
prolonged decay or rise of the count rate; and ( Hi) sources with 
slow modulation or complex variability that could be due to su- 
perimposed flaring events. 

3.1. Flares 

Seven of the TMC sources show evident flares in their light 
curve. Four of them ( V4 1 X7, V8 1 9 Tau, HBC 427 and BP Tau) 
show the typical flare behaviour, characterized by a rapid in- 
crease of the count rate and of the plasma temperature, followed 
by a slower decay. Their light curves and the time evolution of 
the plasma parameters are shown in Figs. Q]|3] The other three 
sources (JH 108, V892 Tau and V410 Tau) show peculiar flar- 
ing variability, i.e. flares that are characterized by a gradual rise 
phase with a flat top and/or a nearly symmetrical shape. Their 
light curves and evolution of the plasma parameters are shown 
in Figs.|4]and|5] 

3.1 .1 . V41 X7 (XEST-23-050/24-042) 

The light curve of V410 X7 (XEST-23-050/24-042, M0.75, 
WTTS) shows two flares that occurred after an initial quiescent 
period lasting for ~ 65 ks. Unfortunately, the XEST-23 exposure 
stopped just after the start of the rise phase of the first flare, while 
at the beginning of the XEST-24 exposure the flare was already 
decaying; therefore the flare peak was not observed. The flare 
is very strong, with an increase of the count rate by more than 
an order of magnitude. The decay phase is initially rapid, with 
an e-folding timescale of ~ 2.1 ks; then a slower decay is ob- 
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HBC 427 = XEST-26-072 - Bin = 1200 s 
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Fig. 2. Left: Light curve and spectral fitting results for HBC 427 (XEST-26-072). Panels are the same as in Fig.Q] Note that in this 
case we show the combined MOS1+MOS2 light curve, the PN being not available. Right: Evolution of the best-fit temperature and 
EM of the hot component for HBC 427. To better clarify the trend with time of the parameters, the points are connected by dotted 
lines. The thick solid line shows the linear fitting to points #3 to #7, corresponding to the flare decay; for clarity they are marked 
with open squares. 



served, after a small bump occurring about 5 ks after the start of 
the XEST-24 exposure. A second smaller flare, with a decay e- 
folding time of ~ 20 ks, is superimposed on the late decay phase. 
As mentioned in Sect. I2.2I our choice of the block size does not 
allow us to derive the plasma characteristics of this second flare 
separately from the decay of the first one. 

The quiescent emission has temperatures of ~ 9 and 24 MK 
with similar EMs. The first flare reached a temperature > 60 MK 
at the peak, as estimated from the best-fit temperature of block 
#2 where the flare is already decaying. Taking this value as peak 
temperature and the decay time of 2.1 ks, and applying to this 
flare the formula for a freely-decaying loop (Eq. (Q}), we can 
put a lower limit on the loop semi-length of L > 5 x 10 10 cm, 
corresponding to ~ 0.4/?+. 

3.1.2. V819Tau (XEST-23-074/24-061) 

V819 Tau (XEST-23-074/24-061) is a binary WTTS with spec- 
tral type K7. It showed a small impulsive flare at the end of the 
XEST-23 exposure, with a factor of 2 increase in the count rate, 
whose decay continues in the XEST-24 observation. The decay 
phase has an e-folding time of ~ 10 ks. The flare is not very hot, 
with T2 = 20 MK in block #2, not much higher than the quies- 
cent value of ~ 15 MK, although the true peak temperature is 
likely hotter since this time interval includes also a significant 
part of the decay phase. We also note that the cool component 
remained constant, within the errors, for the whole observation. 
We can obtain a rough estimate of the loop semi-length using 
the temperatures and EMs of the hot component in blocks #2 
and #3, from which we derive £ ~ 0.9 and L ~ 7 x 10 10 cm, 
corresponding to 0.5 R*. 



3.1.3. HBC 427 (XEST-26-072) 

During the XMM-Newton observation, the K7 WTTS HBC 427 
(XEST-26-072) underwent a strong long-duration flare that cov- 
ered nearly the entire exposure time of ~ 130 ks. After a short 
initial quiescent phase, the count rate increased by a factor of 
5 in ~ 6 ks, then it decreased slowly, returning to a level simi- 
lar to the preflare one after ~ 70 ks. The decay e-folding time 
is 20.2 ±1.0 ks. The time evolution of the plasma parameters 
follows the typical behaviour observed in flares, with the tem- 
perature peaking during the rise phase and the emission measure 
peaking at the flare peak. The quiescent emission has tempera- 
tures of ~ 8 and 25 MK with equal EMs. While the temperature 
and EM of the coolest component does not change significantly 
during the flare, the hot component reaches a temperature as high 
as ~ 72 MK; therefore, for the flare analysis we will assume that 
the flaring plasma is described by the hot component. 

In the right panel of Fig. [2] we show the evolution of the 
flare in the T vs y/EM plane. The maximum loop tempera- 
ture, derived from Eq. (0, is ~ 168 MK. From the linear fit- 
ting of the points from block #3 to block #7 we find a slope 
f = 0.67 + 0.09. This value indicates the presence of signif- 
icant residual heating after the initial ignition. Combining the 
observed decay time with the above value of from Eq. (|2|i we 
derive L = 2.4 + 0.4 x 10 11 cm, i.e. ~ 2 stellar radii. 

From the EM of the hot component in block #2 we can es- 
timate the mean electron density in the loop at the peak of the 
flare, assuming a loop with typical aspect ratio R\ oop /L = 0.1, as 
commonly observed in solar flares: we find n e ~ 2 x 10 10 cm" 3 , 
consistent with the values found for the COUP flares with simi- 
lar loop size dFavata et al.ll2005l) . 
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BP Tau = XEST-28-100 - Bin = 1 200 s 
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Fig.3. Left: Light curve and spectral fitting results for BP Tau (XEST-28-100). Panels are the same as in Fig.Q] Right: Evolution of 
the best-fit temperature and EM of the hot component. See Fig.[2]for details. For the sake of clarity, we plot only the points from #2 
to #9. The slope has been fitted between points #4 and #8. 



3.1.4. BP Tau (XEST-28-100) 

BP Tau (XEST-28-100) is a K7 CTTS. Its light curve shows a 
flare beginning ~ 15 ks after the start of the observation, with an 
increase of the count rate by a factor of ~ 2.5 in ~ 5 ks and a total 
duration of ~ 35 ks. The e-folding decay time is 10.4 + 0.7 ks. 
After the flare the X-ray emission shows low-level variability un- 
til the end of the observation. As for HBC 427, BP Tau shows the 
typical flare behaviour of the plasma parameters, with the tem- 
perature peaking in the rise phase before the EM. The quiescent 
emission before the flare has temperatures of ~4 and 18 MK 
with EM2/EM1 =1.5; the temperature of the cooler component 
increases to ~ 9 MK during the flare, while the hotter compo- 
nent reaches ~ 50 MK. Both EMs increase significantly during 
the flare. 

The right panel of Fig. [3] shows the evolution of the flare 
temperature and EM. From Eq. we derive a maximum loop 
temperature of ~ 100 MK. The slope between points #4 and #8 is 
£ = 0.83 ±0.13, indicating that residual heating is present during 
the decay. Using Eq. (f2]i we derive L = 1.2 + 0.1 X 10 11 cm, 
comparable to the stellar radius. The estimated density at the 
flare peak is « e ~ 3 x 10 I() cirT 3 , similar to the value found for 
HBC 427. 



3.1.5. JH108 (XEST-1 7-066) 

JH 108 (XEST-17-066, WTTS, spectral type Ml) shows a sym- 
metrical flare, with a gradual increase of the count rate by a fac- 
tor of ~ 4 in ~ 8 ks, followed by a decay with similar duration; 
the decay e-folding time is 3.3 ks. This source has quiescent tem- 
peratures of 8 and 30 MK and EM ratio of ~ 1.5. As before, 
the cool component remains nearly constant, while for the hot- 
ter component the temperature is higher during the rise phase, 



reaching 45 MK, and the EM peaks in the second time block, 
corresponding to the peak and decay phases. A rough estimate 
using Eq. gives LS6x 10 10 cm, i.e. ~ 0.7 iJ*. 

3.1.6. V892 Tau (XEST-23-047/24-040) 

V892 Tau = Elias 1 (XEST-23-047/24-040) is a triple system 
composed of a Herbig Ae star, a lo w-mass companion o f spec- 
tral type ~M2 at 4.1" (Elias 1 NE, iLeinert et all 1 1997b and a 
recently-dis covered close com panion at 0.05" with a mass of 
1.5-2 Mr* (ISmith et al.ll2005h . 

For the first 95 ks of the observation this source showed low- 
level variability, with temperatures in the range T\ =5-10 MK 
and T2 = 20 - 35 MK. Then the source underwent a strong flare, 
with a gradual rise lasting for ~ 7 ks and ending in a plateau, at a 
level a factor of ~ 10 higher than the pre-flare count rate, where 
the emission remained nearly constant for at least 10 ks, until 
the end of the observation. For the spectral analysis during the 
flare, the parameters of the cool component were kept fixed since 
they were not constrained, and the flare evolution was entirely 
dominated by the hot component. The plasma heating is strong 
in the first block of the rising phase, where the temperature in- 
creases from ~ 25 to ~ 65 MK; then the temperature stays at a 
nearly constant level, within the errors, around ~ 95 - 100 MK in 
the following three blocks, and starts to decrease only when the 
count rate reaches the plateu in the last time block. The emis- 
sion measure on the other hand continues to increase reaching 
its maximum in the last block. The observed detailed evolution 
of the temperature and EM is a clear example of the flare heat- 
ing process, in which, following the energy release, accelerated 
electrons precipitating into the chromosphere rapidly heat the 
plasma, and produce an evaporation of material which gradually 
fills the loop increasing the emission measure. 
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The long rise time suggests that a large coronal structure is 
involved in the flare. This flare has been studied previously by 
iGiardino et al.l d2004l) . who modeled the rise phase using de- 
tailed hydrodynamical simulations, finding a loop size of the 
order of 1 x 10 cm, equal to ~ 0.5 iJ*, and that a magnetic 
field in excess of 500 G is required to confine the flaring plasma. 
Based on the position of the X-ray source and on the comparison 
with a Chan dra observation where the primary and Elias 1 NE 
are resolved, Giardino et al. (2004) concluded that the flare oc- 
curred in the corona of the Herbig Ae star, which would imply 
the existence of a convective zone in the outer layers of the star, 
generally not expected in stars of this spectral type. However, 
the discovery of a lower-mass close companion suggests that the 
observed X-ray emission and the flare more likely originate from 
it rather than the Herbig star itself. This is supported by the high 
temperatures (20-3 MK) observed in this source and in other 
Herbig Ae/Be stars dStelzer et al.ll2006bl) . similar to those com- 
monly found for later-type PMS stars but significantly higher 
than those (- 5-6 MK) found in the Herbig stars HP 163296 
and AB Aur (ISwartz et al.ll2005b iTelleschi et alj|2006al) . 

3.1.7. V410Tau (XEST-23-032/24-028) 

The WTTS triple system V410 Tau (XEST-23-032/24-028) 
shows a very peculiar light curve with a nearly symmetrical tri- 
angular shape. The count rate rises gradually during the first 
55 ks of observation until block #8, where it reaches a peak a 
factor of ~ 3 higher than at the start of the observation. Then 
the emission decays until the last block; a second, weaker flare 
is superimposed on the decay ~ 9 ks before the end of the obser- 
vation. Assuming that the first 18 ks of the observation represent 
the quiescent level, we find that the rising phase has an e-folding 
time of ~ 20 ± 1 ks. After the peak the emission decays ini- 
tially with an e-folding time of 8.1 + 0.4 ks until the end of the 



XEST-23 exposure; at the beginning of the XEST-24 exposure 
the emission stays at a steady level for 23 ks (block #13), then 
the decay restarts, and continues until the end of the observation, 
interrupted only by the small flare. 

During the observation, the temperature of the cool compo- 
nent does not vary significantly, and the flare evolution is deter- 
mined by the variations of and of the EMs. In the rise phase 
the hot temperature initially increases from 20 to 30 MK, reach- 
ing a maximum in block #4, then decreases again in the follow- 
ing two time intervals; a sudden increase to ~45 MK occurs in 
block #7, just before the peak. The EMs of both plasma compo- 
nents increase during the entire rise phase, reaching their maxi- 
mum at the flare peak, as commonly observed in solar and stellar 
flares. After the peak the temperature and the EMs decrease re- 
turning to the quiescent value. A new temperature increase is 
observed in block #15, corresponding to the second small flare. 
The long rise time, combined with the slow change of the plasma 
temperature during the rise phase, might indicate that we are ob- 
serving a rotationally modulated flare, or the superposition of a 
flaring event and an underlying modulated emission. 

We also note that the best-fit colum density increases by a 
factor of 4 during the flare, following a trend similar to that of 
the EMs. However, we found equally acceptable fits, although 
with slightly higher values (by at most 10%) of the reduced 
X 2 , by keeping Nh fixed to 1 x 10 20 crrT 2 , as found in the first 
time intervals, without significant changes in the time evolution 
of the other parameters. Since this star is a WTTS, which is 
not expected to possess circumstellar material that could justify 
changes in the intervening absorption, and that such low values 
of A^h are at the sensitivity limit of EPIC, we are not confident 
on the significance of the observed column density variations. 

The right panel of Fig.[5]shows the evolution of the temper- 
ature and EM of the hot component during the flare. The slope 
during the decay, between blocks #8 and #14, is f = 0.47 ± 0.12, 
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Fig. 5. Left: Light curve and spectral fitting results for V410 Tau (XEST-23-032/24-028). Panels are the same as in Fig.Q] Right: 
Evolution of the best-fit temperature and EM of the hot component. See Fig.|2]for details. For the sake of clarity, we plot only the 
points from #7 to #15. The slope has been fitted between points #8 and #14. 



indicating that the decay of the flaring structure is largely deter- 
mined by the heating decay timescale. In this case the loop size 
is poorly constrained, given the large uncertainty on we obtain 
in fact L ~ 4 + 3 x 10 10 cm, or 0.1 - 0.4/?*. 

3.2. Very long decays or increases of the count rate 

A group of six sources in our sample showed a monotonic de- 
crease or increase of the count rate over the entire observation. 
One of them, XZ Tau, has a rising light curve, while for the other 
five sources (V830 Tau, HQ Tau, FS Tau, DH Tau and L155 1 5 1) 
a long-lasting decline of the count rate was observed. The light 
curves and best-fit parameters of this group of sources are shown 
in Figs.[6]and[7] 

3.2.1. XZ Tau (XEST-22-047) 

XZ Tau (XEST-22-047) is a binary star composed of two CTTS 
of spectral type M2 and M3.5. This source shows a long mono- 
tonic rise of the count rate by a factor of ~ 4 in ~ 40 ks. An anal- 
ysis of this observa t ion ha s been reported bv lFavata et alj d2003l) 
and lGiardino et all (I2006I) . 

The quiescent emission at the beginning of the observation is 
rather hot, with temperatures of 8 and 40 MK and EM 2 /EMi ~ 
0.7. The cool component does not vary significantly, except for 
an increase of the EM in the last block. The hotter temperature 
peaks at ~ 70 MK in blocks #2 and 3, although the increase is not 
much significant due to the large errors, while the EM continues 
to increase as the count rate increases up to the last block. Our 
results are in agreement with those reported by iGiardino et al.l 
d2006l) . We note that the trend of T 2 and EM 2 resembles that 
observed during the first part of the rising phase of V410 Tau, 
suggesting that we are observing a similar event. 



3.2.2. V830 Tau (XEST-04-016) 

V830 Tau (XEST-04-016, WTTS, spectral type K7) shows a 
long-lasting decay by a factor of 3 in ~ 30 ks. The two temper- 
atures remain nearly constant, within the errors, around ~ 7 and 
20 MK during the entire decay (apart from a possible increase 
in block #3), and the observed decrease of the count rate is due 
only to the decrease of the two EMs with time. We note that also 
the average plasma temperature (weighted with the EMs) does 
not vary during the decay. It is possible that we are observing 
the final stage of a long-lasting flare, when the temperature has 
already returned to the pre-flare level while the EM is still chang- 
ing. An alternative explanation might be rotational modulation: 
the relative amplitude of the count rate variation around the av- 
erage count rate is ~ 60%, compatible with the range of 20- 
70% found in the Orion Nebula Cluster dFlaccomio et al.ll2005l) . 
However, given the high amplitude of the variation, and the fact 
that the 30 ks exposure covers only one te nth of the star's rota - 
tional period (P I0t = 2.75 d, see Table 10 in lGudel et aT]|2006al) . 
we believe that this interpretation is unlikely, although we cannot 
draw definitive conclusions from the available data. 



3.2.3. HQ Tau (XEST-09-026) 

HQ Tau (XEST-09-026) is a close binary WTTS. Its light curve 
shows a long-lasting decay, with a decrease of the count rate by 
a factor of 2 in 30 ks. The spectral analysis shows a significant 
decrease of the hotter temperature, from 34 to 20 MK, and of 
both EMs, while the temperature of the cool component remains 
constant. This behaviour suggests that we are observing the end 
of the decay of a long-duration flare. The e-folding timescale of 
the observed decay is ~ 45 ks. From the decay of the temperature 
and EM of the hottest component we find a high slope ( ~ 2, 
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curve and spectral fitting results for TMC sources showing a long-lasting rise (XZ Tau = XEST-22-047) or decay 
XEST-04-016, HQ Tau = XEST-09-026, FS Tau = XEST-1 1-057). Panels are the same as in Fig.Q] 



compatible with a freely-decaying loop. Using Eq. ([T) we obtain 
L ^ 7 x 10 1 1 cm. Since the stellar radius of this star is not known, 
we cannot determine whether this represents a compact or a large 
loop. 



3.2.4. FS Tau (XEST-1 1-057) 

FS Tau (XEST-1 1-057) is a binary composed of two CTTS of 
spectral type M0 and M3.5. Its light curve decreases by a fac- 
tor of ~ 6 during the first 20 ks, with an e-folding timescale of 
15 ± 2 ks, and then remains at a low quiescent level for the re- 
maining 20 ks. In this case the spectra are well described by just 
one temperature, the low-temperature component being uncon- 
strained due to the high absorption (A^h ~ 1 — 1.6 X 10 22 cm -2 ). 
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Fig. 7. Same as Fig.|]for DH Tau (XEST-15-040) and L1551 51 (XEST-22-089). 



The temperature does not change significantly, decreasing from 
40 MK to 35 MK from the beginning to the end of the observa- 
tion. On the other hand the EM decreases significantly; a small 
but significant decrease of the absor ption is also observ ed. In this 
case, it is not possible to apply the iReale et al.1 (119971) method, 
since the rate of decrease of the temperature is too low (£ = 0.2). 
Applying the formula for a freely-decaying loop we estimate an 
upper limit to the loop length of L ~ 2.7x10" cm(~ 4^ + ).Note 
that, although the peak temperature of the flare is likely higher 
and therefore the above upper limit is underestimated, we expect 
the true loop size to be much smaller because of the presence of 
strong sustained heating, as indicated by the shallow slope (. 

3.2.5. DH Tau (XEST-15-040) 

DH Tau (XEST-15-040) is a binary CTTS with a separation of 
2.3", unresolved by XMM-Newton. It shows a monotonic decay 
by a factor of ~ 3 over the 30 ks exposure. As for V830 Tau, 
the two temperatures are steady during the decay, at ~ 9 and 
~25 MK, while the two EMs decrease. In this case, contrary 
to V830 Tau, also the ratio EMo/EMx decreases, implying that, 
as time proceeds, less and less hot material contributes to the 
emission. This suggests that we are observing the end of a long- 
lasting flaring event. As for the previous source, very strong sus- 
tained heat i ng is present (£ = 0.2), preventing the use of the 
IReale et al] (1 19971) method, and we estimate L < 4.5 x 10 11 cm 
~3.5fl*. 



3.2.6. L1551 51 (XEST-22-089) 

L1551 51 (XEST-22-089, WTTS) also shows a decay by a factor 
of 2 during the first 30 ks, then stays at a constant level for the 
remaining 20 ks of observation. Also for this source the temper- 
atures do not vary, and the observed variation is due only to the 
changing EMs. As for DH Tau, also the EM ratio decreases, im- 



plying a decrease of the average plasma temperature, likely due 
to a flare decay. This source also shows very strong sustained 
heating = 0.16). From the decay e-folding time of 9.5 ks we 
estimate L S 1 X 10 11 cm, comparable to the stellar radius. 

3.3. Other types of variability 

The remaining three XEST sources identified with TMC mem- 
bers show significant variability in the form of slow modulation 
(DN Tau and HD 283572, Fig. [8j or complex flare-like variabil- 
ity (Hubble 4, Fig. . 

3.3.1. DN Tau (XEST-1 2-040) 

DN Tau (XEST- 12-040, CTTS, spectral type M0) shows a mod- 
ulated light curve with a higher count rate at the beginning and 
the end of the exposure and a minimum in the second time block; 
the amplitude of the modulation with respect to the mean level 
is ~ 15%. The plasma parameters do not show significant varia- 
tions, being consistent within the errors, although there is a ten- 
dency for higher Ti and EM\ in the first and last block. 

3.3.2. HD 283572 (XEST-21-039) 

HD 283572 (XEST-21-039) is a G5 WTTS with P mt = 1.55 d; 
the present observation, with an exposure time of ~ 45 ks, there- 
fore covers about one third of the rotational period. The light 
curve shows a clear modulation, with two characteristic levels of 
emission. During the first 20 ks the count rate increases slightly, 
then in block #4 it decreases reaching a level a factor of ~ 1.3 
lower, that is maintained until the end of the observation. The 
two temperatures do not vary significantly, with the cooler one 
staying around 8-9 MK, and the hotter one varying between 
20 and 25 MK with a trend that closely follows that of the light 
curve. On the other hand, the two EMs show significant varia- 
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Fig. 8. Light curve and spectral fitting results for the TMC sources DN Tau (XEST-12-040) and HD 283572 (XEST-21-039), showing 
a slow modulation. Panels are the same as in Fig.Q] 
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tions, in particular in the first part of the observation, where the 
cool EM decreases and the hot one increases, leading to an in- 
crease in their ratio and therefore in the relative contribution of 
the hotter plasma to the emission. These characteristics suggest 



that the observed variability is likely due to the appearing and 
disappearing, due to the stellar rotation, of active regions con- 
taining hotter plasma with respect to the rest of the corona. 



3.3.3. Hubble 4 (XEST-23-056/24-047) 

The light curve of Hubble 4 (XEST-23-056/24-047, WTTS, 
spectral type K7) shows significant low-level flare-like variabil- 
ity superimposed on a slow modulation. The average count rate 
shows a very slow decrease during the first 50 ks, then it in- 
creases by a factor of ~ 1 .5 at the beginning of block #4, and 
decays again slowly until the end of the observation. The two 
temperatures do not vary significantly, having average values of 
~ 9 MK and ~ 20 MK, with the exception of the last block where 
the hot component is slightly higher. As in other sources, the 
variations of the two EMs are more significant, showing a mod- 
ulation with time of both values, but keeping their ratio nearly 
constant, between 1 and 1.2. The observed slow variation of the 
mean count rate might be due to rotational modulation of the 
emission from active regions rotating in and out of view. 



4. Sources not associated with known TMC 
members 

In addition to the TMC members, we have studied three ad- 
ditional XEST sources (XEST-05-031 = HD 283810, XEST- 
16-031 = 2MASS J04195676+2714488, and XEST-22-024 = 
HD 285845) that are likely non-members of the TMC but that 
are strong X-ray sources and show significant variability. Their 
light curves and spectral parameters are shown in Figs. [TOl and 
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4.1. HD 283810 (XEST-05-031) 

XEST-05-03 1 is identified with HD 283810, a K5V star with Ha 
emission, which is probably an older foreground star since it has 
radial velocity inconsistent with TMC membersh ip and does not 
show significant Li absorption (Herbi g"et al.lll98a) : assuming a 
main-sequence object, its photometry (V = 10.74, B — V = 1.03) 
locates it at a distance of ~ 60 pc. Its light curve shows a flare at 
the end of the observation, with an increase of the count rate by a 
factor of 2.5 and an e-folding decay time of ~ 2.4 ks. The hotter 
temperature rises from ~ 9 MK during quiescence to 18 MK in 
block #3 and reaches a maximum value of 26 MK in block #4, 
after the flare peak. Both EMs increase significantly at the flare 
peak. Given the irregular trend of th e temperature and EM, we 
cannot apply the iReale et al] d 19971) method; using the formula 
for a freely-decaying loop and the observed temperature in block 
#3, we derive a loop semi-length of L < 3 x 10 10 cm, which 
is comparable t o the stellar radius for a K5 main-sequence star 
dSiess et al . 200(J. 

4.2. HD 285845 (XEST-22-024) 

XEST-22-024 is identified with HD 285845, which is a fore- 
ground binary system at a distance of 90 pc; the absence of the 
Li absorption line at 6708 A indicates that the system is not com- 
posed of PMS stars dWalter et alJll988tlFavata et al.ll2003l) . The 
mean properties of the X-ray spectrum have been discussed by 
Fava ta etail (|2003), who reported a spectral fitting with a 2-T 
model characterized by high Ca and Ne abundances with respect 
to a Fe abundance around 0.26 Z . 

The light curve of this source shows a complex variability. A 
flare occurred at the beginning of the exposure, with a gradual 
rise lasting ~ 5 ks, followed by a decay on which several minor 
impulses are superimposed. The initial decay just after the peak 



has an e-folding time of ~ 2 ks. The emission reaches a steady 
level ~ 30 ks after the peak. As observed in other sources, the 
cool temperature does not vary significantly during the flare evo- 
lution. The hotter temperature peaks in block #2, during the rise 
phase, at 40 MK, and remains steady at a level of ~ 30- 35 MK 
up to block #7, when it decreases to ~ 15 — 25 MK, remain- 
ing at this level until the end of the observation. The observed 
light curve variations are mostly due to the EMs of both compo- 
nents, that vary significantly throughout the observation. Given 
the complexity of the light curve, it is not possible to derive the 
loop size by fitting the T vs EM decay. We can make a rough 
estimate assuming the initial decay time of 2 ks and Eq. ([1), ob- 
taining L ~ 3 x 10 10 cm. A ssuming a G6 main-sequence star, 
this corresponds to ~ 0.4 tf* dSiess et alj|2000l) . 

4.3. 2MASS J041 95676+271 4488 (XEST-1 6-031) 

Source XEST-16-031 has an IR counterpart in the 2MASS cata- 
logue, 2MASS J04195676+2714488, located at ~ 0.2" and with 
magnitudes J = 12.38, H = 11.80 and K = 11.54 mag; an- 
other fainter 2MASS source is present at 4.6", however the X- 
ray source is most likely associated with th e former one, given 
the accuracy of t he XEST positions (see iGiidel etail [2006a; 
IScelsi et al.lS006h . 2MASS J04195676+2714488 has photom- 
etry inconsistent with TMC membership, being located on the 
main sequence in the colour-magnit ude diagrams for the dis- 
tance of the TMC ( Scel si et ai1l2006l) . and is therefore likely to 
be a field late-type star. The uncertainties in the photometry al- 
low a distance between 80 and 190 pc for a main sequence star. 

During the first 23 ks, the source was in a quiescent state with 
very low emission. Unfortunately in this part of the observation 
(blocks #1 to 3) the background was very high and dominating 
the observed count rate (see Fig. [TTJ, preventing the possibility 
of performing spectral analysis. In block #4, ~ 24 ks after the 
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Fig. 11. Left: Light curve and spectral fitting results for XEST-16-031. Panels are the same as in Fig.Q] Right: Evolution of the 
best-fit temperature and EM of the hot component. See Fig.[2]for details. The slope has been fitted between points #10 and #14. The 
slope during the decay of the first small flare (points #4 and #5) is also indicated. 



start of the observation, a flare occurred, followed by a second 
much stronger flare 7 ks later. The first flare has a peak count 
rate a factor of ~ 5 higher than the quiescent level, while the 
second flare increases the count rate by more than one order of 
magnitude. The two flares have very similar decay times, with 
e-folding times of ~ 1.6 and 1.5 ks, respectively. This suggests 
that they may have occurred in loops of comparable size. 

The absorption is very low, consistent with a column density 
A^h = 1 X 10 20 crrT 2 , and to better constrain the fit parameters it 
was kept fixed to this value. The spectrum was well fitted with 
only one temperature component. Since the emission during the 
quiescent phase is negligible with respect to the flare phase, we 
assume that the 1-T model is representing essentially the flaring 
plasma during the time blocks from #4 to #14. 

The right panel of Fig. [TT] shows the evolution of the tem- 
perature and the EM during the observation. In block #4, cor- 
responding to the peak of the first, small flare, the plasma tem- 
perature is 22 MK; the temperature and EM decrease in block 
#5 during the flare decay, and remains steady in the following 
time interval. In block #7, at the start of the rise of the second 
flare, the temperature increases to ~ 42 MK; a second re-heating 
is observed in block #9, corresponding to the flare peak, where 
the EM reaches its maximum value. The flare cooling phase from 
point #10 to point #14 proceeds with a slope £ = 0.42 ±0.07, that 
indicates the presence of significant sustained heating after the 
initial ignition. It is worth noting that the slope between points 
#4 and #5, relative to the decay of the first small flare, has a very 
similar value of 0.37. Using Eq. (O we derive L ~ 1 x 10 9 cm 
for the first flare, and L ~ 5 ± 4 x 10 9 cm for the second flare. 
These values are significantly smaller than those estimated for 
the other sources. Although the precise nature of this source can- 
not be assessed here, given the very low absorption we suggest 
that it could be an older M-type star located just in front of the 



cloud and showing flare events analogous to those observed on 
the Sun and active late-type stars. 

5. Discussion and conclusions 

In this paper we have studied a sample of 19 bright variable X- 
ray sources detected in the XEST survey. Our sample includes 16 
known TMC members (10 WTTS, 5 CTTS and a Herbig Ae star) 
plus three additional sources unrelated to the cloud but show- 
ing flaring events. The studied sources show different types of 
variability, in the form of flares, either with fast rise and slower 
decay or with symmetrical shapes, continuous rise or decays, 
slow modulation possibly due to rotation, and complex variabil- 
ity, with flare-like events superimposed on a slow modulation 
or decay. Using detailed time-resolved spectroscopy we have in- 
vestigated the changes of the parameters of the emitting plasma 
(temperatures and emission measures) and of the hydrogen col- 
umn density, and, in the case of flares, we have derived infor- 
mations on the size of the involved coronal structures. The main 
parameters derived for the studied sources are summarized in 
Tabled 

The quiescent emission has typical temperatures Ti ~ 4 - 
10 MK and T 2 ~ 15 - 35 MK, with EM 2 jEM x ~ 0.7 - 2, 
consistent with the values found in other studies of young PMS 
stars (e.g lFeigelson & Montmerielll999l: IWolk et alj|2005l) . No 
significant difference is evident in our small sample between 
the spectral characteristics of CTTS and WTTS: both the ranges 
of the parameters and their median values are similar for the 
two classes. The Herbig Ae star V892 Tau has spectral char- 
acteristics very similar to those found for TTS stars, in agree- 
ment with other studies dHamaguchi et alj 120051 : IStelzer et alj 
2006b), supporting the suggestion that the emission might come 
from its cool, unresolved close companion rather than from 
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Table 3. Summary of the parameters derived for the XEST sources. For all sources we give the temperature of the hot component 
(T q ) and the luminosity (Lxq) in the quiescent (or lowest count rate) intervals; a range is given for sources with slow modulation or 
complex variability. For sources showing flares or prolonged rise/decays, we also give the peak temperature (r pea k) and luminosity 
(LxpeakX the rise and/or decay e-folding timescales (T r i se and t^c), and the derived loop size L. 



XEST ID 


Optical ID 


T 


Lxq 




^Xpcak 


^"rise 


T dec 


L 


L/R* 






(MK) 


(10 30 ergs-') 


(MK) 


(10 30 ergs-') 


(ks) 


(ks) 


(10" cm) 




04-016 


V830 Tau 


22 


2.8 


>33 


>8.5 




11.4 






09-026 


HQTau 


19 


4.6 


>34 


>9.1 




45.2 


>7.1 




11-057 


FS Tau 


36 


1.4 


>42 


>6.3 




15.2 


<2.7 


4.1 


12-040 


DN Tau 


19-27 


1.0-1.4 














15-040 


DH Tau 


20 


4.7 


> 24 


> 9.8 




34.2 


S4.5 


3.5 


17-066 


JH 108 


28 


1.5 


46 


3.7 


3.8 


3.3 


£0.6 


0.7 


21-039 


HD 283572 


19-26 


10-14 














22-047 


XZTau 


40 


0.9 


70 


3.8 


35.3 








22-089 


L1551 51 


16 


1.4 


>17 


>2.5 




9.5 


<1.1 


1.1 


23-032/24-028 


V410 Tau 


20 


3.8 


43 


11.7 


20.9 


8.1 


0.1-0.7 


0.1-0.4 


23-047/24-040 


V892 Tau 


23 


9.5 


100 


46.8 


3.4 




1.0" 


0.5 


23-050/24-042 


V410 X7 


24 


0.7 


>63 


>4.9 


1.5 


2.1 


>0.5 


0.4 


23-056/24-047 


Hubble 4 


19-30 


3.5-5.4 














23-074/24-061 


V819 Tau 


14 


2.0 


19 


3.8 


1.3 


10.1 


0.7 


0.5 


26-072 


HBC 427 


25 


3.2 


72 


15.8 


3.2 


20.2 


2.4 


1.9 


28-100 


BP Tau 


18 


1.5 


47 


4.9 


2.6 


10.4 


1.2 


0.9 


05-031 


HD 283810 


10 


1.5* 


26 


3.2 6 


0.5 


2.4 


0.3 




16-031 


J04195676+27 






21 


3.5 


0.6 


1.6 


0.01 












42 


20.0 


0.8 


1.5 


0.01-0.09 




22-024 


HD 285845 


18 


2.1^ 


40 


6.5 c 











": from Giardino et al. (2004) 
: computed for d = 60 pc 
' : computed for d = 90 pc 



the Herbig Ae star itself. We find that in most sources the 
cool plasma component does not vary significantly, even dur- 
ing flares, while the observed time evolution can be totally as- 
cribed to variations in the hot component. Similar results have 
been obt ained for the O rion PMS stars studied in the COUP 
survey (|Wolk et aill2005l) . as well as for older active stars (e.g. 
lAudardetal.1120011) . 

Nine of the studied sources (including two TMC non- 
members) show evident flaring activity. The flares observed on 
TMC members have total duration between ~20 and 70 ks 
(~ 5-20 hr), with e-folding rise timescales of ~ 1 - 4 ks, except 
for the peculiar flare on V410 Tau with ~ 20 ks, and decay 
timescales of ~ 2 - 20 ks. An additional 5 sources show gradual 
decays over 30-50 ks, that might represent the decay of long- 
lasting flares, as suggested by the decreasing plasma temperature 
and emission measure, and one source shows a prolonged rise 
(with e-folding timescale of ~ 35 ks) with spectral characteris- 
tics similar to those observed during the long rise phase of the 
atypical flare on V410 Tau. Unfortunately, the typical 30-40 ks 
exposure times of the XEST survey introduce a bias against the 
detection of long-duration flares in our observations. Indeed, we 
detected flares with total duration of ~ 50 - 70 ks in the archival 
fields having exposure times of ~ 100 ks. It is conceivable that 
the sources with gradual decays were undergoing flares of sim- 
ilar duration or even longer, such as found in the C OUP survey, 
where events lasting up to 3 days were observed dFavata et alj 
l200HlWolk et~aT]l2005l) . 

Spectral analysis of the flaring sources shows peak temper- 
atures from 40 MK up to 100 MK for the strongest flares, and 
peak luminosities between 4 x 10 30 and 5 x 10 31 erg s _1 . On the 
other hand, the flares observed on the two TMC non-members 
show significantly shorter rise and decay times ($ 1 ks and 



~2 ks, respectively), and peak temperature in the range 20- 
40 MK. 

For four of the flaring sources (the WTTS HBC 427 
and V410 Tau, the CTTS BP Tau, and the non-member 
2MASS J04195676+2714488) we could perform a detailed 
analys is of the decay phase using the method by iReale et all 
dl997l) : the method was applied also to V819 Tau, for which 
however only two points in the decay are available. In all these 
cases we find that significant residual heating must be present 
during the decay, governing the observed light curve evolution. 
For the other flaring sources, we do not have enough intervals 
in the decay phase to apply the method, and we could only 
estimate upper limits to the loop size using the formula for a 
freely-decaying loop. The loop size is fully constrained only for 
the CTTS BP Tau, which has L = 1.2 x 10 11 cm, comparable 
to the stellar radius, and for the WTTS HBC 427, which has 
L = 2.4 x 10 11 cm corresponding to ~ 2/?*; the latter star 
has the longest decay time, and hence the longest loop size, 
among the sources for which the entire flare evolution is ob- 
served. The other TMC members have loop lengths in the range 
4 - 7 x 10 10 cm, smaller than or comparable to the stellar radius. 
For the sources with continuous decay, only one (HQ Tau) has 
a steep slope, compatible with a freely-decaying loop with no 
additional heating after the initial ignition; for this star we find 
L ~ 7 x 10 11 cm. The other four sources have upper limits be- 
tween 1 - 5 x 10 11 cm, comparable to the stellar radius for the 
two WTTS, but equal to ~ 4fl* for the two CTTS (FS Tau and 
DH Tau). Finally, for the TMC non-members we find loop sizes 
in the range 1 x 10 9 - 3 x 10 10 cm. 

We mention that flare characteristics similar to those found 
here for TMC members are observed a lso for the CTTS SU Aur 
(XEST-26-067. lFranciosini et alJl2007l) . that showed three flares 
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during the observation with rise and decay times of ~6 ks 
and ~ 5 - 9 ks, respectively, and peak temperatures of ~ 50 - 
140 MK. The data do not allow a detailed flare ana l ysis, h owever 
using the parameters reported by Franciosin i et al.l (120071) we de- 
rive upper limits to the loop semi-length of ~ 1.1 — 1.8x 10 11 cm, 
comparable to the stellar radius. 

Previous observations of PMS stars have shown the presence 
of compact flaring structures, w ith L £ R+, similar to what is ob- 
served in active late- type stars dFavata et alJ|200lUGrosso et all 
2004tlGiardino et al.l2006tlGudel et alj2004tlReale et al.ll2 004). 



Favat a et al. (2005) analyzed a sample of strong flares detected 
in the COUP survey, finding both compact structures, of size 
shorter than a stellar radius, and very extended structures, with 
lenghts of 5 - 207?+. Such long structures possibly represent 
magnetic loops connecting the stellar surface with the circum- 
stellar disk. Our sample shows generally loops of size compara- 
ble to or smaller than the stellar radius: the longest loop with a 
fully constrained size is about 2 stellar radii in length, which is 
compatible with a loop anchored on the stellar surface. We stress 
that this loop has been observed on a WTTS, that should not pos- 
sess a circumstellar disk, and that a similar loop size (~ 1.6 R+) 
has be en found also on the WTTS V827 Tau by iGiardino et al.l 
(2006). A possible hint for large loops of size ~ 4i?* is found 
for the two CTTS with long-lasting decays, suggesting that these 
stars might have indeed loops connecting the star and the disk; 
however the estimated size is highly uncertain, since we are ob- 
serving only a small part of the decay, therefore we cannot draw 
any definitive conclusion on the size of the emitting structures. 

Finally, we have studied two sources showing possible ro- 
tational modulation with amplitudes of ~ 15%, and two other 
sources (one unrelated to the TMC) with flare-like variability su- 
perimposed on a slow modulation or on a flare decay. Except for 
the last case, we do not find significant variations in the plasma 
temperatures, and the observed variability is mainly determined 
by variations of the emission measures. 
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Table 2. Best fit values of the spectral model parameters. Errors are 90% confidence ranges for one interesting parameter; where errors are not 
given, the parameter was held fixed to the tabulated value. F x and L x are the unabsorbed X-ray flux and luminosity in the 0.3-7.3 keV band. L x is 
computed from F x assuming a distance of 140 pc for all stars, except for HD 283810 and HD 285845, where we used the photometric distance of 
60 and 90 pc, respectively. Note that for J04 195676+27 (XEST-16-031) the reported luminosity must be considered only indicative, since the true 
distance of this star is not known. 
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30.52 






2 


9 19+0.38 
z - -0.29 


0.86^1? 


9 14+O.57 
Z - -0.31 


0.78+»;« 


2.63+oj 2 


1.12 


120 


1.57 


30.57 






3 


9 qo+0.46 

Z -"^-0.38 


n 7 a+0.08 
-0.08 


9 43+1.36 
z -^-'-0.62 


0.59t! 2 


0.89!»;» 9 


1.23 


97 


0.64 


30.18 



21-039 HD 283572 



1 


0*6™ 


75 +0 - 03 
-0.03 


, 04+0.22 

-0.20 


5 48+ 129 

J>lo -1.06 


6 37+0- 63 

"■ J '-0.68 


0.99 


347 


5.46 


31.11 


2 


0.72+0.05 


74+0-01 
"■ '^-0.02 


2 03+ - 07 
z - -0.10 


4 61+0- 40 
HD -0.66 


9^+0.38 
S - zj -0.24 


1.05 


825 


6.05 


31.15 


3 


59+ 04 
UJ -0.02 


76 +0 01 
"■' -0.01 


9 99+O.O8 
z - zz -0.12 


3 34+0.33 
J - J +-0.45 


8 96 +0 35 


1.09 


782 


6.21 


31.16 


4 


0.64+0;! 8 


70+ 05 

"■'"-0.06 


1 oo+0.19 
l- 8i -0.16 


3.86+| : « 


'■ J "-0.66 


0.99 


220 


5.39 


31.10 


5 


0.68+o-oa 


75 +0 - 02 
"■' J -o.oi 


1 O7+0.07 
l- S '-0.08 


4 14+0- 56 
H- -0.38 


6 23 +0 - 26 

°- ZJ -0.28 


1.05 


715 


4.91 


31.06 


6 


64 +0 06 
"■" -0.06 


73+0.02 

-0.03 


1 67+0- 06 

1D -0.08 


3 48+ 047 


6 11+040 


0.98 


616 


4.50 


31.02 


7 


66 +0 07 

"■""-0.07 


74+0- 02 
"■' -0.02 


1 o«:+0.09 
180 -0.08 


4 12+0- 59 
^• 1Z -0.66 


5 O7+0.38 
J -^'-0.31 


1.16 


553 


4.74 


31.05 



1 


9 45+O.6I 


n 7^+0.08 

"■' J -o.n 


3 44+2.31 
-0.95 


45+O.I2 
"■ HJ -0.09 


0.31+0;0 9 


1.44 


75 


0.36 


29.93 


2 


2 05 +0 - 47 


0.82+0-15 


5 79+4.49 
J -' z -1.71 


49+0 16 
"■ H7 -0.12 


0.68+0-J 2 


0.90 


80 


0.65 


30.18 


3 


, qi+0.39 


0.77!0.;3 


^ nq+2-06 

6.03_ 151 


45 +o.l4 


1 95+O.I6 
■ J -0.09 


1.04 


143 


1.00 


30.37 


4 


9 99+O.23 
Z ' ZZ -0.21 


0.78+0-" 


4 49+O.64 
-0.45 


0.49+O.J 2 


1 £-3+0.12 


0.95 


278 


1.30 


30.49 


5 


2 56 +0 - 29 

—0,26 


0.79:0-09 


3 50 +0 - 52 
■ -0.42 


0.82+o.J 9 


9 90+0.21 
Z - ZJ -0.16 


0.94 


230 


1.60 


30.58 



22-089 


L1551 51 


1 


0.87+0-21 


65+ 014 

"■" J -0.17 


1 50+ 049 

1 - J "-0.23 


1 01+ 031 
1 -" 1 -0.26 


1 53+O.24 
■ -0.47 


0.92 


63 


1.08 


30.40 






2 


1.00+0|? 


62+ 06 

"■"-0.15 


1 25 +0 - 59 
1 - ZJ -0.21 


1-13^6 


0.61+0; 2 i 


1.19 


90 


0.72 


30.23 






3 


0.62*$ 


63+ 05 
"■"-'-0.06 


1 -3,7+0.63 
• J ' -0.23 


0.85+0.J2 


0.38+0;! 4 


0.95 


101 


0.59 


30.14 
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YTJCT TT~\ 
ADJ 1-1D 


uptical ID 


rSlOCK 




T 
1 1 






ElVl2 


„2 


dot 










u 
it 


iu cm 


Kev 


Kev 




cm 3 






-2 -1 

erg cm s 


erg s"i 


23-032+24-028 


V410 Tau 


l 


0.10 


74+ 002 

-0 02 


1 70+0- 14 

-0 08 


1 7Q+0.12 
-0 14 


2.00+0|4 


1.27 


685 


1.63 


30.58 






2 




o.75-;; : :;! 


2 05+° : 25 
-0.15 


1 90+0 .14 

-0.16 


2.80+019 


1.07 


540 


2.10 


30.69 






3 


°- 09 -006 


0.75+004 


2 30+0-22 

^- JU _0 26 


1 90+0 .16 

-0 12 


oq+0.31 
Jm -0 24 


1.18 


498 


2.43 


30.76 






4 


0.15^ 


0.79+0-03 


2.48+o:23 


9 C9+O+9 
-0 19 


3 60+0 26 


1.06 


494 


2.84 


30.82 






5 


0.19*$ 


0.77+004 


2.33+0-g 


3 01+031 

-0.28 


4 04 +042 


1.08 


373 


3.25 


30.88 






6 




0.75+003 


2 - 10 t?2 


o co+0.31 
-0.24 


' u 42 


1.02 


571 


3.63 


30.93 






7 




0.82+0-^ 


3.67+1-66 


4 47 +0.63 


4.42+o17 


1.00 


144 


4.41 


31.02 






8 




0.77+0-04 


2.50+og 


n ot+0.42 
j.oj_q ^2 


6.70+059 


0.87 


355 


4.95 


31.07 






9 




0.78+0-04 


2 -°2!o+6 


o 9<r+0.40 
■ J - z -'_0 38 


5.67+0-52 


1.04 


369 


3.99 


30.97 






10 


0.37*»£ 


0.77!004 


2 07+0-43 
-0 21 


3.27+g ; 4g 


4.35+049 


0.88 


292 


3.40 


30.90 






11 


21+ 007 
"■-0.06 


76+ 03 

"•'"-0.04 


1 90+02 1 

1 - 7 "-0.17 


3 03+0- 28 

J - -0.28 


1 99+0.35 
J - ZZ -0.35 


1.15 


377 


4.27 


31.00 






12 


^ ™ ^_i_n no 

0.20+0$ 


u -' -0.04 


1 QQ+0.34 

-0.23 


2.66+03J 


9 '7Q+0.38 
z - -0.40 


1.00 


291 


2.38 


30.75 






13 


o.iC 

— u.UJ 


73 +0 - 02 

— U.U1 


1.79+ 011 

—0.10 


!- 98+ 009 

— o.oy 


2.52+0-07 


1.15 


818 


1.94 


30.66 






14 


0.07+006 


0.77+0-03 


j 75+0.I6 


1 55+0- 16 

— — 28 


2.09+o:35 


1.03 


448 


1.58 


30.57 






15 


0.12 +0 |° 


0.79+005 


2^9+0.54 


1.74+0 28 


2.66+035 


1.03 


212 


2.06 


30.69 






16 


0.10 


0.58+0-|2 

— U. 1 / 


i ro+0+8 

n on 
— U.zU 


0.99+JI 


2.78+0-38 

— V.jj 


1.19 


93 


1.58 


30.57 






17 


0.10 


0.76+0;04 


1.81+0-32 
— U. lo 


1.41+0 19 

—0.16 




0.91 


311 


1.41 


30.52 


23-047+24-040 


V892 Tau 


1 


8.96+1; 67 


0.82+0-22 


2.98+o»4 


3.15+278 


a QQ+0.63 
J,j:7 _0.73 


1.10 


254 


3.13 


30.87 






2 


10.51+144 


0.69+0-3? 


2-931^3 


5.86+2,4 


4 QQ+1.25 

-1 67 


0.79 


77 


5.08 


31.08 






3 


8.90+i ; io 

-0.77 


0.92+0-;; 

-0.16 


2.87+0-52 

-0.36 


3 25+1 27 

-1.08 


4.00+075 

-0.73 


0.97 


338 


3.46 


30.91 






4 


8.14+^ 


1.06+0-1* 


3.02+0-™ 


2.89+|-86 


4 QQ+103 

-1 18 


1.05 


194 


3.86 


30.96 






5 


9.59 + [;20 


0.81+0-J5 


9 q-7+0.40 

-0 27 


4.59! 2 -54 


4 19+0:73 

-0.80 


1.13 


242 


3.99 


30.97 






6 


1 9 1-66 


0.38+o : o 7 


1 84+° :i8 
1 - -0.16 


13.90+i2 8 30 


c 1 t+0.66 
-0.59 


0.97 


232 


7.52 


31.25 






7 


1105 185 


"■ -0.08 


1 Q^+ l '-^' z - 

1-7 -0.20 


6.00_353 


3.69+0;6i 


1.05 


295 


4.09 


30.98 






8 


10 - 73+ 067 


0.49 


5.71+2-04 


6.00 


6 -21 + 061 


1.03 


90 


6.38 


31.18 






9 


12.07+J-jg 


0.49 


8.04+406 


6.00 


13.71+j5j 


1.02 


85 


11.60 


31.44 






10 




0.49 


8.62+412 


6.00 


18.04+] : 48 


0.72 


124 


14.59 


31.54 






11 


11 53+ 079 

■ -0.78 


0.49 


8.23+2-99 


6.00 


21.38+148 


0.99 


181 


16.72 


31.59 






12 




0.49 


7.24+0-65 


6.00 


26-72!°;73 


0.99 


862 


20.08 


31.67 


23-050+24-042 


V410 X7 


1 


7.41+f^ 


0.77+0-3' 


2.07+2-38 


0.331J;60 


0.35+009 


1.16 


130 


0.30 


29.85 






2 


8.19+i$ 


1 36+0- 68 

-0.36 


5.43+2-36 


2.42+(-58 

— 1,10 


1.69+J-76 

— l.UO 


0.92 


105 


2.07 


30.69 






3 


8.34+(4' 

— 1.0/ 


0.68+0-28 


2.55+049 

— U.4/ 


0.63+1-51 


2.07+0-56 

—U.J 1 


0.97 


139 


1.29 


30.48 






4 


10.46 +4 -8i 

— J .8 / 


0.68+0;24 


1-78!144 


2.02+209 

— 1.51 


1.34+1 69 

—0.73 


1.06 


116 


1.45 


30.53 


23-056+24-047 


Hubble 4 


1 


9 -39+0.09 
J -0.09 


74+0 02 

"' -0.02 


1.81+0|2 


9 C4+O.I6 
Z --' -0.16 


9 49+0.24 
z -+^-0.19 


1.19 


680 


2.18 


30.71 






2 


1 ^Q+0 37 
1- -032 


77 +, '- u/ 
"■"-0.10 


~ rin _|_n 1(\ 

2.28+0;76 


i CT+0 38 
1.5j„ 35 


1 x 4-0 40 
1 -'"-0.38 


0.96 


89 


1.50 


30.55 






3 


9 i r+0.20 
z - -0.19 


76+0-04 
"•'"-0.05 


1-84!^ 


2 07+0-3 1 
z '"'-0.26 


2-07!»;» 


1.00 


239 


1.80 


30.63 






4 


9 94+0.13 
Z - Z# -0.13 


77+0-03 

"■"-0.04 


2.08+0-28 


2 30+0-26 


2.80+024 


1.00 


443 


2.28 


30.73 






5 


Z - Uy -0.16 


76+0-04 

"■'"-0.04 


1.87+0-22 


9 19+0.26 
z - lz -0.24 


9 1 4+0-28 
^- -0.28 


1.06 


302 


1.85 


30.64 






6 


9 is+0.14 
z - -0.13 


n 7Q+0.04 
"• -0.12 


1 63+0-36 


1.86+031 


1 88+0-35 
i.oo_ 033 


1.16 


494 


1.60 


30.58 






7 


9 9Q+0.44 
z -^-0.28 


73+0-05 
■ -0.10 


2 63+ 1 - 15 

z -° -0.58 


9 47+0-42 
z - -0.38 


1 7Q+0.42 
1,/7 -0.38 


0.92 


165 


1.94 


30.66 



23-074+24-061 V819Tau 


1 


2 11+0 23 
z - -0.20 


40+0-02 

"■+"-0.02 


1 31+0-13 

^ -0.07 


1 41+O.26 
1 -^ 1 -0.19 


92+0 .07 

"■ ?z -0.07 


1.06 


438 


0.93 


30.34 




2 


9 94+ -58 
-0.51 


n 30 +0.09 


lu "-0.19 


1.60+085 


9 qo+0.28 
2- 3S - .35 


0.87 


80 


1.63 


30.58 




3 


, 79+0.49 
!-' z -0.42 


a 37+O.O6 
"•-"-0.05 


1 99+O.I2 

!- ZZ -0.11 


1.08+049 


1 1 O+0.14 
L1S -0.16 


0.89 


100 


0.89 


30.32 




4 


2.23+046 


36+0-03 

"■ -0.04 


1 20+0-12 

1 - z "-0.12 


1 34+0.47 
1 - -0.31 


89+0 .12 

"■ oy -0.12 


0.93 


218 


0.87 


30.31 



18 E. Franciosini et al.: X-ray bright variable sources in the Taurus Molecular Cloud 

Table 2. continued 



1-1D Uptical ID 


rSlOCK 




T 
1 1 




LM\ 


T7A/f 
LIV12 


„2 


aot 


17 /in - 12 






# 
ft 


iu cm 


Kev 


Kev 


in'3 

1U 


cm 3 






-2 -1 

erg cm s 


ergs -1 


26-072 HBC 427 


l 


20+ - 20 

-0 19 


a 7 q+0.06 
-0 07 


2 14+0-90 

43 


1 51+0- 26 


1 48+ - 35 
■ -0 33 


0.88 


104 


1.34 


30.50 




2 


29+ 032 

"■ -0.28 


93+ 016 
"■ -0.18 


6-18^1 


2 00+ 92 
z - -0.82 


5 30+0^94 

-0.82 


0.99 


73 


4.32 


31.01 




3 


°- 57+ a23 


j-, — r. _i_n i "7 


3.31+0;^ 


1.81+0; 92 


10.89+ ^ 4 


0.95 


132 


6.70 


31.20 




4 


°- 49+ 026 
— U.zO 


0.89+0JI 


2 - 88+ om 

— U.04 


2.68+J- 36 


6.68+! :3 « 

— l.oz 


0.94 


74 


4.62 


31.04 




5 


— U.lo 


—U.UO 


2.43+0- 3 5 

— u.zy 


2.16+0-40 

— U.jJS 


S.lO+og 


1.01 


190 


3.41 


30.90 




6 


°- 51 -012 


°- 74+ 005 

—U.LO 


2.01+0-38 
— U.zU 


2.23+0- 3 ' 


3.27+0|| 

— U.43 


1.31 


220 


2.45 


30.76 




7 


0.28+0- 

—U.lo 


0.83+015 

—U.UO 


1.89+0-80 

— U.43 


2.23+0-4' 

— U.43 


1.83+0- 63 


1.10 


128 


1.82 


30.63 




8 


032 -a!o 


°- 74+ 006 


1 - 54+ 14 


1.76 +0 f 




1.15 


217 


1.57 


30.57 


28-100 BPTau 


1 


°- 9 ^a22 


a -17+0.04 
-0 03 


1 57+0- 11 

-0 11 


0.63+0|4 


0.94+005 


1.25 


226 


0.63 


30.17 




2 


1.00 


41+° :u 

-0 06 




0.68tj4 


1 9^+0.16 
J -0.19 


0.76 


75 


0.87 


30.31 




3 


1.00 


79+ a24 
-0 23 


4.04+|-'5 


0.78+0- 2 6 


1 oa+033 


0.98 


79 


1.46 


30.54 




4 


1 26+ 012 


75+° :,)6 

-0 06 


2 65 +°- 23 

Z - UJ _0 17 


0.89+0|4 


q 9Q+0.21 
-0 21 


0.93 


419 


2.06 


30.69 




5 


1 30+ 029 


n 70+O.O6 

-0 09 


2 00+049 


1 55+0- 28 


1 93+0^52 


1.23 


133 


1.59 


30.57 




6 


O.SS+i" 4 


n 77+O.O8 
-0.09 


9 it+0.47 
• -0.30 


1.06*& 


1 55+ 031 


0.80 


101 


1.17 


30.44 




7 


0.95! 2 ;« 


n 70+O.O6 
U - /8 -0.08 


i o o +0.24 
1,iS2 -0.19 


0.80!0|6 


1 97+O.I9 
i,z,/ -0.21 


1.14 


156 


0.90 


30.33 




8 


0.64_ 038 


n -70+0 21 
U - /8 -0.08 


1 ^n+0 57 

1.60 Q15 


n co+0 16 
°- 52 -0.19 


n n/1+0 16 
°- 94 -0.26 


0.97 


103 


0.62 


30.16 




9 


i 97+0.26 
-0.20 


34+0-02 
""^-0.02 


1 50+ 007 
1 - J "-0.08 


66+° 14 


92+0 05 

"■'^-0.05 


0.99 


345 


0.63 


30.17 




10 


1.00 


40+0- 05 


1 62+0- 25 

1 - DZ -0.12 


54+0 05 
"■ J ^-0.07 


94+0.07 

"■^-0.09 


1.20 


179 


0.60 


30.15 




11 


0.87+oj?, 


37+ -03 

"■ J '-0.02 


i +0.09 
^ -0.07 


40+0.07 
U ' 4y -0.07 


71+0 05 

"■ 11 -0.05 


1.05 


348 


0.48 


30.05 




12 


°- 69 -02, 


n 38+0.05 


1 O7+0.21 
L87 -0.22 


49+0.14 

"•^'-0.09 


1 01+0 .07 
1 -" 1 -0.07 


1.03 


218 


0.64 


30.18 




13 


1 03+ - 35 

1 - u -'-0.40 


31+0- 07 

-0.03 


I 99+0.17 
1 " 4y -0.07 


45+0 .19 
"■ -0.16 


94+0.07 

"•^-0.21 


0.98 


198 


0.56 


30.12 


05-031 HD 283810 


1 


°- 01 -ooi 


0.41+0-05 


0.98+008 


2.16+0- 4 5 


1 95+O.28 

—0 28 


1.21 


165 


1.66 


30.22 




2 


O-OO-mo 


0.41+0-02 


0.79+0-05 


1-69+0^ 


1 93+031 

26 


1.49 


420 


1.48 


30.17 




3 


0.01 


0.60+0* 

— U.UO 


!- 56 -o?i 


3.01+0- 4 5 

~u.4y 


4.66+059 

— U.D / 


1.07 


172 


3.24 


30.51 




4 


0.01 


0.69+006 


2 - 27+ 072 


■3 CT+0.54 
78 


L72+ 68 


1.18 


88 


2.32 


30.37 




5 


0.01 


0.19+0-0' 

— U.Uo 


0.75+0-05 

— U. 1 J 


0.96+' 4 :f 


4.00+ 2 '' 

— U.4 / 


1.32 


106 


1.89 


30.28 


16-031 J04 195676+27 


4 


11+ 014 
w - 11 -0.11 




1. 85+0-g 




3 34+0.21 

-0.21 


1.10 


110 


1.48 


30.54 




5 


0.10 




1 65+0- 17 

-0 11 




1 si+0.09 
1 -° 1 -0.07 


1.29 


91 


0.78 


30.26 




6 


0.10 




1.63+^ 




9 3^+0.19 

-Q02 


1.50 


123 


1.74 


30.61 




7 


0.10 








6.56.^8 


0.87 


69 


3.71 


30.94 




8 


0.10 




9 qo+0.39 




9.88!0|? 


1.02 


68 


4.80 


31.05 




9 


0.10 




3.02+0-30 




16 20+ 1 - 15 


1.12 


161 


8.58 


31.30 




10 


0.10 




2.01+O-g 




13.33+0;68 


1.05 


65 


6.08 


31.16 




11 


0.10 




1.92+0-g 




11.22+o ; 56 


1.17 


78 


5.05 


31.07 




12 


0.10 




1.63+0- 22 




7 - 06 -009 


1.51 


102 


3.02 


30.85 




13 


0.10 




1.64+0" 

— U. 1U 




3.81+0 [« 

—U.lo 


1.26 


100 


1.63 


30.58 




14 


0.10 




1.24+006 

—U.UO 




1 ro+0.07 
n no 

— u.uy 


1.47 


95 


0.65 


30.18 


22-024 HD 285845 


1 


°- 42 -o!o 


64+ 003 


9 <=9 +0.54 
z " -0.30 


q 97+O.3I 
-0 24 


q 97+O.33 

-0 38 


0.93 


295 


3.00 


30.66 




2 




(\ 77+O.O6 
_0 06 


•5 aq+6.30 

-1 28 


4.09+0^0 


a a 9 + 1.01 
- , -- ,z, _0 87 


0.96 


86 


3.58 


30.74 




3 


0.47+006 


061 +o:o4 

-0 16 


2.86+0^ 


3 50+ 42 

—0 38 


S 9Q+0.63 
—0 47 


1.09 


232 


4.21 


30.81 




4 


0.36+«i» 


66 +0 06 


3 10+ 1 - 43 


a aq+0.54 


3 50+0 68 


1.01 


119 


3.30 


30.70 




5 


°- 40 -o!i 




3.01+^ 


3.20+0^ 


2 - 70+ 033 


1.09 


215 


2.78 


30.63 




6 


39+° 09 


65+ o: 05 


3.04+0-48 


2.99+ , 28 


3.62+0- 33 


0.90 


339 


3.19 


30.69 




7 


51+ 009 

"■ J -0.08 


n A/1+0.03 


i £,+0.33 


2 70+0- 40 
z - /u -0.24 


a C7+0.31 

J - J '-0.45 


0.91 


368 


2.67 


30.61 




8 


34 

"■ -0.22 


,, ,_ + n no 
fi7 

""'-0.05 


2 - 48 -o:s6 


tit 45 
3 1 3+ u + ; ' 
J - -0.42 


1 t * +0 49 
'• -0.59 


0.71 


105 


1.99 


30.48 




9 


32+°- 12 
UJ -0.12 


n fio+o.03 
0.68 003 


i C9+0.32 


9 70+O.24 
^■' J -0.24 


1 11+° 24 
1 ' 11 -0.26 


0.95 


191 


1.64 


30.40 




10 


20+ - 20 
"■^"-0.20 


66+ 06 
"■""-0.58 


i ^7+0.39 
■ -1.09 


2.02+0- 3 5 


1 11 + 1 - 39 
1-11 -0.31 


0.67 


65 


1.32 


30.30 






59+ 09 
U ' J '-0.I1 


66+ 02 

"■""-0.06 


| O9+0.29 
1 -^ z -0.24 




2 16+0- 28 


0.94 


217 


2.16 


30.52 




12 


°- 38 -oli 


71+ 03 

"•"-0.04 


1 84+0-65 
1 -° -0.37 


2 61+0- 28 

Z - Oi -0.35 


1 36+0- 40 

1 -' O -0.35 


1.04 


230 


1.72 


30.42 




13 


35+ 05 
-0.06 


66+0° l 

"■""-0.02 


1 56+0- 14 
i.jo_ 13 


9 to+0.09 
JO -0.09 


89+o .09 

"■ -0.09 


1.10 


428 


1.38 


30.32 



E. Franciosini et al.: X-ray bright variable sources in the Taurus Molecular Cloud 



List of Objects 



'V410X7' on page|4 
'V819 Tau' on pagef' 
'HBC 427' on pageg 
'BP Tau' on pagejn 
'JH 108' on page@ 
'V892 Tau' on page[l 
'V410 Tau' on page[' 
'XZ Tau' on page|| 
'V830 Tau' on page[l 
'HQ Tau' on pageR 
'FS Tau' on pagef 
'DH Tau' on pagellOl 
'L1551 51' on page[TC 
'DN Tau' on pageflC 
'HD 283572' on page 
'Hubble 4' on pagefll 
'HD 283810' on page 
'HD 285845' on page[ 
'2MASS J04195676+2714488' onpage[l2 



